Pentaquark baryon production at the Relativistic Heavy Ion Collider 
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Production of pentaquark baryons in central relativistic heavy ion collisions is studied in a 
kinetic model. Assuming that a quark-gluon plasma is produced in the collisions, we first determine 
the number of O"*" produced from the quark-gluon plasma using a parton coalescence model, and 
then take into consideration its production and absorption in subsequent hadronic matter via the 
reactions KN ^ 6, KN ^ irO, and nN ^ KQ. We find that although the final number is 
affected by hadronic interactions, it remains sensitive to the initial number of produced from 
the quark-gluon plasma, particularly in the case of a small width as imposed by the N and 
K'^d scattering data. Because of small baryon chemical potential in the hot dense matter produced 
in these collisions, the number of produced anti-0 is only slightly smaller than that of . 

PACS numbers: 25.75.-q,25.75.Dw,25.75.Nq 



I. INTRODUCTION 



In recent experiments on nuclear reactions induced by 
photons ,1,, iSj ! kaons Q , and neutrinos 3] j pro- 

duction of baryons consisting of five quarks uudds has 
been inferred from the invariant mass spectrum of K'^n 
or K'^p. The extracted mass is around 1.536 GeV with a 
width in the range 20-25 MeV except Ref.f^ which gives 
a width of 9 MeV. All these widths reflect the resolu- 
tion in the experiments, so the actual width of 0+ is 
expected to be smaller. The observed properties of this 
pentaquark baryon are consistent with those of the Q'^ 
baryon with spin J ^ 1/2, isospin / = 0, and strangeness 
S" = -f 1 that was originally predicted by the chiral soli- 
ton mod el i8 i| and recently studied using the Skyrme 
model IglllolT ll l lT ^ ll ^ lTil. the constituent quark model 
[ll|l3jl3Uhe chiral quark model T?,']^, the QCD sum 
ruleslialil |22J, and the lattice QCD 2^ J4j. Stud- 
ies of the Q'^ production mechanism in these reactions 
have also been carried out |2^ |2^ |23, . Since both 
kaon and nucleon numbers are not insignificant in the 
hadronic matter formed in relativistic heavy ion colli- 
sions, the O"*" may also be produced in these collisions. 
Using a statistical model, which assumes that the 0+ 
is in chemical equilibrium with other hadrons, Randrup 
[2^ has estimated its abundance and finds that there is 
about one per unit rapidity in central Au-f Au colli- 
sions at y^SNN — 200 GeV available from the Relativistic 
Heavy Ion Collider (RHIC). As the quark-gluon plasma 
is expected to be formed in the initial stage of heavy ion 
collisions at RHIC, and it was suggested that formation 
of the quark-gluon plasma would enhance the production 
of hadrons consisting of strange quarks [3(ij | , it is of inter- 
est to know if this is also the case for 0+ production and 



to determine its production from the quark-gluon plasma 
relative to that from the hadronic matter. 

In this Letter, 0"*" production in central heavy ion col- 
lisions at RHIC is studied in a kinetic model that starts 
from the final stage of the quark-gluon plasma, goes 
through a mixed phase of quark-gluon and hadronic mat- 
ters, and finally undergoes hadronic expansion. The pro- 
duction of 0+ from the quark-gluon plasma is modeled by 
the coalescence model, which has been shown to describe 
quite well not only the particle yields and their ratios 3l| 
but also their transverse momentum spectra ,32, 33, 3^. 
Production and absorption of O"*" in the hadronic mat- 
ter are taken into account via the reactions KN <-^- 0, 
KN ^ ttQ, and ttN ^ KQ. We find that the number 
of produced from the quark-gluon plasma is appre- 
ciable and the final number after the hadronic phase 
remains sensitive to the initial number of produced 
from the quark-gluon plasma. Furthermore, a slightly 
smaller number of its antiparticle G)~ is also expected to 
be produced at RHIC as a result of small baryon chemi- 
cal potential in the matter produced from these collisions. 
We note that the quark coalescence model has also been 
used recently to study in heavy ion collisions at RHIC the 
elliptic flow of 0+ i.e., the azimuthal anisotropy of 
their momentum distribution in the plane perpendicular 
to the beam directions. 



II. COLLISION DYNAMICS AT RHIC 

To model the dynamics of central relativistic heavy ion 
collisions after the end of the quark-gluon plasma phase, 
we use the boost invariant picture of Bjorken ,36j aug- 
mented with accelerated transverse expansion. Specifi- 
cally, the volume of produced fireball is taken to evolve 
with the proper time according to j37l | 
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V{t) =tt[Rc+ vc{t - Tc) + a/2(T - tc)^ 



(1) 



where Rc — 8 fm and tq — 5 fm/c are final transverse 
and longitudinal sizes of the quark-gluon plasma, corre- 
sponding to a volume of 1,006 fm'"* in central Au-fAu 
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collisions at y^i^vw = 200 GeV. Taking the critical tem- 
perature to be Tc = 175 MeV and a transverse flow veloc- 
ity vc = 0.4c, the total transverse energy of quarks and 
gluons in the midrapidity {\y\ < 0.5) is then about 1,067 
GeV, if we take quarks and gluons to be massive with 
TOg = 500 MeV, rriu = mj. ^ 300 MeV and = 475 
MeV in order to take into account the nonperturbative 
effects of QCD near critical temperature 38], and to have 
chemical potentials fib = 10 MeV and fig = to ac- 
count for observed final antibaryon to baryon ratio at 
RHIC. The above total transverse energy includes a bag 
energy of about 133.8 GeV based on a bag constant of 
about 133 MeV/fm'^, which is determined from the pres- 
sure difference between the quark-gluon plasma and the 
hadronic matter at Tc- Requiring that final hadronic 
matter freezes out at temperature Tp = 125 MeV and 
has a transverse fiow velocity of 0.65c, as extracted from 
measured hadron spectra in these collisions, then leads 
to a freeze out volume V{tf) « 11, 322 fm^^ if we assume 
that the fireball expands isentropically. The resulting to- 
tal transverse energy of final hadronic matter is 788 GeV 
and is comparable to that of midrapidity hadrons mea- 
sured in experiments |39f . It is, however, less than the 
initial transverse energy of midrapidity quarks and glu- 
ons. This is consistent with the reduction of midrapidity 
hadrons and their total transverse energy as a result of 
hadronic rescatterings seen in studies based on the trans- 
port model ^3 • Because of the small pressure near phase 
transition and in hadronic matter (411] . acceleration in 
the transverse expansion is chosen to have a small value 
a = 0.02 c^/fm in order to obtain a lifetime of the ex- 
panding matter comparable to that from the transport 
model. 

As the quark-gluon plasma expands, it gradually con- 
verts to hadrons leading to a mixed phase of quark-gluon 
plasma and hadronic matter at constant temperature Tc- 
The fraction of hadronic matter /h (t) can be determined 
by requiring that the total entropy Stot of the fireball re- 
mains constant during the phase transition [i^, i.e., 

/H(r).H(Tc) + (1 - /H(r)).QGP(rc) = (2) 

V(t) 

where sh and sqgp are the entropy densities of hadronic 
matter and quark-gluon plasma, respectively, which we 
treat as noninteracting ideal gases with finite transverse 
flow. The mixed phase ends at time th ~ 7.5 fm/c 
when /h(th) reaches one, and the resulting pure hadronic 
matter is also assumed to expand isentropically until 
Tp « 17.3 fm/c when its temperature drops to Tp = 125 
MeV. In panels (a) and (b) of Fig^ we show the resulting 
time evolution of the volume and temperature of mid- 
rapidity particles in Au-|-Au collisions at ^/snn — 200 
GeV. 

For normal hadrons such as pions, kaons, and nucle- 
ons, we take them to be in chemical equilibrium with 
baryon chemical potential [Ib = 30 MeV, charge chem- 
ical potential hq = MeV, and strangeness chemical 
potential us = 10 MeV, similar to those of the quark- 
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FIG. 1: (Color online) Time evolution of the volume (a), tem- 
perature (b), and abundance (c) of mid-rapidity particles in 
Au-|-Au collisions at ^snn = 200 GeV. 



gluon plasma. Since the chemical potentials vary weakly 
with the temperature of an isentropically expanding mat- 
ter in heavy ion collisions at RHIC energies we ne- 
glect their time dependence. Evaluating the densities of 
these hadrons using the relativistic Boltzmann distribu- 
tions with inclusion of transverse flow, the time evolution 
of their abundance is shown in panel (c) of FigQ] We 
note that the final numbers of pions, kaons, and nucle- 
ons at freeze out including contributions from decays of 
resonances are comparable to those measured in experi- 
ments. 



III. e+ PRODUCTION FROM THE QGP 

The number of 0+ baryons that are produced from 
the quark-gluon plasma can be estimated using the coa- 
lescence model j32i |. i.e., given by the product of a statisti- 
cal factor gs , which denotes the probability of combining 
uudds quarks into a color neutral, spin 1/2, and isospin 
hadronic state, and the overlap of the quark phase-space 
distribution function fq{xi,pi) with the Wigner distribu- 
tion function /e of 0+. The latter corresponds to the 
probability of converting the above hadronic state into 
and depends on the quark spatial wave functions in 
the O"*". Explicitly, the number is expressed as 

(3) 
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with da denoting an element of a space-like hypersur- 
face. Since the coalescence model can be viewed as for- 
mation of bound states from interacting particles with 
energy mismatch balanced by other particles in the sys- 
tem, neglecting such off-shell effects is reasonable if the 
binding energy is not large as in the case of 0+ and/or 
if the production process is fast compared to the inverse 
of the binding energy. We note that the 8"*" is produced 
from the quark-gluon plasma during its hadronization 
and thus exists only in the resulting hadronic matter. 

To determine the statistical factor ge , we note that the 
quark wave function of 0"*" in the color-spin-isospin space 
can be expressed as a linear combination of all possible 
orthogonal flavor, color, and spin basis states, with coef- 
ficients depending on the quark model used for O^, e. g., 
see Since these coefficients are normalized to one, 

the probability for the five quarks and antiquark uudds 
to form a hadron with quantum numbers corresponding 
to a spin up 0+ is simply given by the probability of 
finding these quarks in any one of these color-spin-isopin 
basis states, i.e., 1/3^ x 1/2^ = 1/7776. Including also the 
possibility of forming a spin-down 0"*" doubles the proba- 
bility. As a result, the probability factor is gs — 1/3888. 

For the phase-space distribution function of quarks, 
we assume that they are uniformly distributed in the 
transverse plane and their momentum distributions are 
relativistic Boltzmannian in the transverse direction but 
uniform in rapidity along the longitudinal direction. Im- 
posing also the Bjorken correlation of equal spatial (77) 
and momentum [y) rapidities, then the quark momentum 
distribution per unit rapidity is 



fqiv,y,P±) 9qSiv - y)exp 



T 



(4) 



where 17^ = 6 is the color-spin degeneracy of a quark. 

The Wigner distribution function of O"*" depends on its 
internal quark wave functions. For simplicity, we take it 
to be Gaussian in space and momentum, i.e., 
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feix;p) ^8^e^p[-J2^-J2 ^l^l ^ (5) 



i=l 



where and are the relative spatial and momentum 
coordinates of the five quarks in and are related to 
Xi and Pi by the normal Jacobian transformation, i.e.. 



= ( 1 + - 



S.7 = l ' 



(6) 



and similarly for the momentum space coordinates. The 
width parameter of for the i-th relative coordinate is 
defined as of — (fiiUj)^^ with 



(7) 



Carrying out the spatial integrations in Eq.lPJ, the 
number of O"^ produced from the coalescence of uudds 
quarks is then 



iVc-: 



\ 5 5 



(47ra2)! 
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-(8) 



In the above, k are the four relative momenta de- 
fined previously but determined in the center-of-mass of 
formed pentaquark, and mi± = ^ mf + are trans- 
verse masses of the quarks. The numbers of u, d, 
and s quarks in the quark-gluon plasma with rapidities 
\y\ < 0.5 are denoted, respectively, by N^, Nd, and Ng; 
and a is the width parameter calculated according to 
the mass niq of the light quark. At re their values are 
Nu = Nd — 245 and Ng ~ 149, if we take into account 
the effect of gluons by converting them to quarks accord- 
ing to the quark flavor composition in the quark-gluon 
plasma as in Ref . '32'| . The width parameter a in the 
0+ Wigner distribution function is related to the size of 
0+. Since the latter is not known empirically, we thus 
choose the value of cr to fit the proton root-mean-square 
radius with the harmonic oscillator wave functions, i.e., 
(7 — (rp)^/2 ~ 0.86 fm 14^, corresponding to a 0+ root- 
mean-square radius 



2\l/2 



<J 1 



4~P^ 



1/2 



0.9 fm. 



(9) 
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Evaluating numerically the transverse mass integrals 
in Eq.® using the Monte Carlo method introduced in 
Ref.,32iJ, we find that the number of 8+ produced in 
one unit of rapidity is about 0.19. Compared to predic- 
tions from the statistical model, this number is smaller 
than that estimated in Ref. by about a factor five and 
that from a more recent analysis by about a factor 
two. The coalescence model also allows us to evaluate 
the number Ny of neutral hyperons A and E'' as well 
as those from decays of S° and S~ that are produced 
from the quark-gluon plasma. The resulting number is 
about 10. Although this is about a factor of two smaller 
than that measured in Au-|-Au collisions at ^/s = 130 
GeV li^l , including contributions from resonance decays 
and production in the hadronic matter is expected to 
bring it closer to the experimental data. In deriving 
Eq.®, we have neglected the effect of transverse flow 
on the quark momentum distributions. Including this ef- 
fect makes the spatial integrals in Eq.lPl more involved. 
Evaluating them by the Monte Carlo method, we find 
that its effect on produced 0+ number is small. It is in- 
teresting to note that the momentum integrals, i.e., the 
last factor in Eq.®, can be evaluated analytically if the 
quark momentum distributions are non-relativistic, lead- 
ing to the expression (1 -I- 2mqTa'^)~'^. This gives a 
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number of 0.64, which is about a factor of 3.4 larger than 
the relativistic case. 



IV. HADRONIC EFFECTS ON 6+ 
PRODUCTION 

The abundance of 0+ can change during subsequent 
evolution of the hadronic matter as a result of hadronic 
absorption and production. In hadronic matter, 0+ can 
be produced as a resonance from interactions of K and N . 
The cross section for this process is given by the Breit- 
Wigner formula. It can also be produced from reactions 
like KN ttQ and TriV KQ. Cross sections for these 
reactions can be estimated by considering KN ttG as 
a w— channel and ttN — > KQ as an s— channel nucleon- 
pole diagram. Using the coupling constant Qknb — 4.4, 
determined from the width of 0+ which we take to be 
Fe = 20 MeV, we have evaluated these cross sections 
with empirical form factors. Details of this calculation 
can be found in Ref.[2^. The 9+ can be destroyed in 
the hadronic matter either by decay, i.e., @ KN, or 
by the inverse reactions ttO ~> KN and KQ ttN with 
cross sections related to those of production cross sections 
via the detailed balance relations. 

We have also evaluated the thermal average of 8"*" pro- 
duction and absorption cross sections. In the tempera- 
ture range 125-175 MeV of interest here, their values for 
production cross sections are about 1 mb for (aKN^ev) , 
0.3 mb for {(Jkn^ttQv) , and 1-3 ^b for (c^Ar^ife^)' 
where v is the relative velocity of two initial hadrons and 
(• • •) denotes the average over their momentum distribu- 
tions. For Q'^ absorption cross sections, they are about 
6 mb for {(Tt^q-^knv) and 0.6 mb for (cxe^TrN^)- The 
cross sections are thus larger for Q^ production reactions 
induced by K than tt and also for Q^ absorption than 
production reactions as a result of large Q'^ mass. 

In terms of thermal averaged cross sections and the 
densities of pions (htt), kaons (ji/f), antikaons («/f), and 
nucleons (jin), the time evolution of Q'^ abundance iVe 
is determined by the kinetic equation 48, 49], 

dNe 



dT 



+ {{aKN->7rev) + {o-KN^ev))nKnNV}i 

-TeNe - {(trq^^nv) n^Ne 

- {(Jttq^knv) n^Ne- (10) 

Since hadronization of the quark-gluon plasma takes a 
finite time of th — tc — 2.5 fm/c, we expect Q^ to be 
produced continuously from the quark-gluon plasma in 
the mixed phase, with a rate proportional to the volume 
of the quark-gluon plasma. This is included in the first 
term i?QGp(T) of above equation. 

Using the time evolution of volume, temperature, and 
hadron abundance shown in Fig^ the resulting time evo- 
lution of 8+ abundance is shown in FigI2Ia). The solid 
curve is obtained with 0.19 initial from the quark- 
gluon plasma as given by the coalescence model discussed 
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FIG. 2; (Color online) Time evolution of O"*" abundance in 
Au+Au collisions at ^/snn ~ 200 GeV for different initial 
numbers of produced from the quark-gluon plasma as well 
as different O"*" width. 



in the above. It is seen that the final 0+ number is en- 
hanced to 0.46. To see how the final 8"^ number depends 
on its number produced from the quark-gluon plasma, we 
have repeated the above calculation with different num- 
bers of initial 8"'" baryons. As shown in Figl^^a), the 
final 8+ number is about 0.41 and 0.62, respectively, 
for (dotted curve) and 0.94 (dashed curve) initial 8+ 
baryons, roughly corresponding to a change of 30% in the 
value of a or the size of 8+ . Our results demonstrate that 
although the final number of 8+ produced in relativistic 
heavy ion collision is affected by hadronic interactions, 
it is nonetheless sensitive to the initial number of 8+ 
produced from the quark-gluon plasma. 

This sensitivity would be even stronger if the width 
of 8+ is smaller than that used here, i.e., Fe = 20 
MeV. In this case, not only the decay probability of 8+ 
is reduced but also the hadronic cross sections become 
smaller as they are proportional to the square of the cou- 
pling constant gKN& , which is proportional to Fe . This 
is shown in Figs. |2Ib), (c), and (d) for the abundance 
of 8+ using different 8+ widths of 10, 5, and 1 MeV, 
respectively. The smaller width of 5 MeV or less is actu- 
ally more consistent with available K^N and K'^d data 
ia HI m m m. with such & smaU 8+ width, the 
hadronic effects become insignificant, and the final 8^" 
abundance is essentially determined by production from 
the quark-gluon plasma. 

On the other hand, the dependence of final 8''" num- 
ber on its initial number from the quark-gluon plasma 
disappears only if hadronic cross sections are large. This 
happens when the latter is increased by about a factor 
of 3, corresponding to a 8+ width of 60 MeV. Although 
such a large 8^ in-medium width is unlikely due to the 
weak coupling of 8^ to KN, the final 8''" number in 
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this case is about 0.35 and is comparable to that from 
the statistical model prediction of pq . 

V. DISCUSSIONS 

In the kinetic equation, Ea. HlU|) . we have not consid- 
ered reactions involving more than three particles such 
as ttQ ^ ttKN. Taking its cross section the same as 
that for the reaction ttQ ^ KN, Ea. (|10|l can be gener- 
alized to include both production and annihilation. 
We find that including this reaction affects the final 8+ 
abundance by less than 10%. Also, medium effects on 
the properties are not included. Because of the weak 
<d~^ coupling to KN and the small hadron density after 
hadronization of the quark-gluon plasma, these effects 
are not expected to be important. 

Our result that the final number of exotic pentaquark 
0+ produced in heavy ion collisions at RHIC is sensi- 
tive to its initial number produced during the hadroniza- 
tion of the quark-gluon plasma is, however, based on 
the assumption that 0+ production from initial nucleon- 
nucleon collisions is not significant compared to that 
from the hot dense matter formed during the collisions. 
This is expected to be the case as string fragmentation 
to multiquark-anti-multiquark pairs is suppressed com- 
pared to fragmentation to diquark-anti-diquark pairs, 
from which normal three-quark baryons are produced in 
high energy proton-proton collisions. Furthermore, both 
the quark flavor structure and mass of 0^, which be- 
longs to the antidecuplet of SU(3), are similar to those of 
the decuplet 17"*" (sss), whose abundance in central rela- 
tivistic heavy ion collisions is known to be more than an 
order of magnitude larger than that expected from initial 
nucleon-nucleon interactions [4^ . 

Our study can be applied to production of anti-0+ in 
heavy ion collisions at RHIC. With the light antiquark 
to quark ratio of 0.89 and antistrange to strange quark 
ratio of 1 as a result of small quark baryon chemical po- 
tential /If, = 10 MeV, which gives an antiproton to proton 
ratio of (0.89)'^ ~ 0.7 at midrapidity consistent with ex- 



perimental observations, the 8~ to 0+ ratio is expected 
to be (0.89)"* ~ 0.63. Relativistic heavy ion collisions at 
RHIC thus also allows us to find the O^, which is not 
likely to be produced in either photo- or kaon-nuclcus 
reactions. 



VI. SUMMARY 

To summarize, we have used a kinetic model to study 
production of pentaquark O"*" baryon in central relativis- 
tic heavy ion collisions by including contributions both 
from the quark-gluon plasma via the parton coalescence 
model and from the hadronic matter through hadronic 
reactions between nucleons and pions or kaons. We find 
that a substantial number of 0+ are produced from the 
quark-gluon plasma. Because of the expected narrow 
width, the O^ interacts weakly in the hadronic mat- 
ter. As a result, the final number of O^ is not much 
affected by interactions in the hadronic matter and is 
thus sensitive to the initial number of 0+ produced from 
the quark-gluon plasma. Study of 0+ production in rel- 
ativistic heavy ion collision thus provides the possibility 
of understanding the dynamics of quark-gluon hadroniza- 
tion. 



VII. ACKNOWLEDGMENTS 

We thank Berndt Miiller and Jorgen Randrup for com- 
ments on an earlier version of this paper. This Letter was 
based on work supported in part by the US National Sci- 
ence Foundation under Grant No. PHY-0098805 and the 
Welch Foundation under Grant No. A-1358. VG was 
further supported by the National Institute of Nuclear 
Physics (INFN) in Italy, while LWC by the National Sci- 
ence Foundation of China under Grant No. 10105008 and 
SHL by the Korean Research Foundation under Grant 
No. KRF-2002-015-CP0074. 



[1] T. Nakano et al., LEPS Collaboration, Phys. Rev. Lett. 

91, 012002 (2003). 
[2] S. Stepanyan et al., CLAS Collaboration, Phys. Rev. 

Lett. 91, 252001 (2003). 
[3] V. Kubarovsky et al., CLAS Collaboration, Phys. Rev. 

Lett. 92, 032001 (2004). 
[4] J. Barth et al., SAPHIR Collaboration, Phys. Lett. B 

572, 127 (2003). 
[5] A. Airapetian et al., HERMES Collaboration, Phys. Lett. 

B 585, 213 (2004). 
[6] V.V. Barmin et al., Phys. Atom. Nucl. 66, 1715 (2003); 

Yad. Fiz. 66, 1763 (2003). 
[7] A.E. Asratyan et al., Phys. Atom. Nucl. 67, 682 (2004); 

Yad. Fiz. 67, 704 (2004). 



[8] D. Diakonov, V. Petrov, and M. Pohakov, Z. Phys. A359, 
305 (1997). 

[9] M. Praszalowicz, Phys. Lett. B 575, 234 (2003). 
[10] M.V. Polyakov and A. Rathke, Eur. Phys. J. A 18, 691 

(2003) . 

[11] H. Walliser and V.B. Kopeliovich, J. Exp. Theor. Phys. 
97, 433 (2003). 

[12] B.K. Jennings and K. Maltman, Phys. Rev. D 69, 094020 

(2004) . 

[13] D. Borisyuk, M. Faber, and A. Kobushkin, 

hep-ph/0307370 
[14] N. Itzhaki, I.R. Klebanov, P. Ouyang, and L. Rastelli, 

Nucl. Phys. B 684, 264 (2004). 
[15] Fl. Stancu and D.O. Riska, Phys. Lett. B 575, 242 (2003). 



6 



[16] M. Karliner and H.J. Lipkin, hep-ph/0307243 

[17] R.L. Jaffe and F. Wilczek, Phys. Rev. Lett. 91, 232003 

(2003) . 

[18] A. Hosaka, Phys. Lett. B 55, 571 (2003). 

[19] L.Y. Glozman, Phys. Lett. B 575, 18 (2003). 

[20] S.L. Zhu, Phys. Rev. Lett. 91, 232002 (2003). 

[21] R.D. Matheus, F.S. Navarra, M. Nielsen, R. Rodrigues 

da Silva, and S.H. Lee, Phys. Lett. B 578, 323 (2004). 
[22] J. Sugiyama, T. Doi, and M. Oka, Phys. Lett. B 581, 167 

(2004) . 

[23] F. Csikor, Z. Fodor, S. D. Katz, and T. G. Kovacs, JHEP 

0311, 070 (2003). 
[24] S. Sasaki, hep-lat/0310014 

[25] W. Liu and CM. Ko, Phys. Rev. C 68, 045203 (2003). 
[26] W. Liu, CM. Ko, and V. Kubarovsky, Phys. Rev. C 69, 

025202 (2004); W. Liu and CM. Ko, ibid., 69, 045204 

(2004); nucl-th/0309023 
[27] Y. Oh, H. Kim, and S.H. Lee, Phys. Rev. D 69, 014009 

(2004). 

[28] S.L Nam, A. Hosaka, and H.C Kim, Phys. Lett. B 579, 
43 (2004). 

[29] J. Randrup, Phys.Rev. C 68, 031903 (2003). 

[30] J. Rafelski and B. Miiller, Phys. Lett. B 101, 111 (1982). 

[31] T.S. Biro, P. Levai, and J. Zimanyi, Phys. Lett. B 347, 

6 (1995); ibid. 472, 243 (2002); Phys. Rev. C 59, 1574 

(1999); T.S. Biro, T. Csorgo, P. Levai, and J. Zimanyi, 

Phys. Lett. B 472, 243 (2002); P. Csizmadia and P. Levai, 

Phys. Rev. C 61, 031903 (2000). 
[32] V. Greco, CM. Ko, and P. Levai, Phys. Rev. Lett. 90, 

202302 (2003); Phys. Rev. C 68, 034904 (2003). 
[33] R.C Hwa and CB. Yang, Phys. Rev. C 67, 034902 

(2003); 064902 (2003). 
[34] R.J. Fries, B. Miiller, C Nonaka, and S.A. Bass, Phys. 

Rev. Lett. 90, 202303 (2003); Phys. Rev. C 68 044902 

(2003). 

[35] C. Nonaka, M. Asakawa, S.A. Bass, R.J. Fries, and B. 
Mueller; nucl-th/03120Mj 



[36] J.D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[37] CM. Ko, X.N. Wang, B. Zhang, X.F. Zhang, Phys. Lett. 

B 444, 237 (1998). 
[38] P. Levai and U. Heinz, Phys. Rev. C 57, 1879 (1998). 
[39] A. Basilevsky et al. (PHENIX Collaboration), Nucl. 

Phys. A715, 486c (2003). 
[40] Z.W. Lin, S. Pal, CM. Ko, B.A. Li, and B. Zhang, Phys. 

Rev. C 64 011902 (2001); B. Zhang, CM. Ko, B.A. Li, 

and Z.W. Lin, ibid. 61, 067901 (2000). 
[41] F. Karsch, Nucl. Phys. A698, 199 (2002). 
[42] K. Kajantie, J. Kapusta, L. McLerran, and A. Mekjian, 

Phys. Rev. D 34, 2746 (1986). 
[43] L. Grandchamp and R. Rapp, Phys. Lett. B 523, 60 

(2001); R. Rapp, Phys. Rev. C 63, 054907 (2001). 
[44] CE. Carlson, C. D. Crone, H.J. Kwee, and V. Nazaryan, 

hep-ph/0312325 
[45] G.G. Simon et al., Nucl. Phys. A333, 381 (1980). 
[46] J. Letessier, G. Torrieri, S. Steinke, and J. Rafelski, Phys. 

Rev. C 68, 061901(R) (2003). 
[47] C. Adler et al. (STAR Collaboration), Phys. Rev. Lett. 

89, 092301 (2002); K. Adcox et al. (PHENIX Collabora- 
tion), ibid. 89, 092302 (2002). 
[48] CM. Ko and L.H. Xia, Phys. Rev.C 38, 179 (1988). 
[49] CM. Ko, V. Koch, Z.W. Lin, K. Redlich, M. Stephanov, 

and X.N. Wang, Phys. Rev. Lett. 86, 5438 (2001). 
[50] R.A. Arndt, I.I Strakovsky, and R.L. Workman, Phys. 

Rev. C 68, 042201(R) (2003). 
[51] J. Haidenbauer and G. Krien, Phys. Rev. C 68, 

052201(R) (2003). ^ 

[52] S. Nussino, hep-ph/0307357 

[53] R.N. Cahn and G.H. Trilling, Phys. Rev. D 69, 011501 
(2004). 

[54] R.A. Arndt, Ya.I. Azimov, M.V. Polyakov, I.I. 
Strakovsky, and R.L. Workman, Phys. Rev. C 69 035208 
(2004). 



